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Abstract 

Ph ' In this short note, we study the local well-posedness of a 3D model for incompressible 

■^^ ■ Navier-Stokes equations with partial viscosity. This model was originally proposed by 

ri , Hou-Lei in [?]. hi a recent paper, we prove that this 3D model with partial viscosity will 

develop a finite time singularity for a class of initial condition using a mixed Dirichlet 

Robin boundary condition. The local well-posedness analysis of this initial boundary 

value problem is more subtle than the corresponding well-posedness analysis using a 

standard boundary condition because the Robin boundary condition we consider is 

►-j>. , non-dissipative. We establish the local well-posedness of this initial boundary value 

CO I problem by designing a Picard iteration in a Banach space and proving the convergence 

^N ■ of the Picard iteration by studying the well-posedness property of the heat equation 

with the same Dirichlet Robin boundary condition. 
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S ■ 1 Introduction 
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In this short note, we prove the local well-posedness of the 3D model with partial viscosity. 
The 3D model with partial viscosity has the following form: 

ut = 2utpz 

ujt = {u'^)z + i^^oj , (x, z) G O = J7x X (0, oo), (1) 

-AV' = oj 

where f^x = (0, a) x (0,a). Let V = {(x, z) | x G ilx; z = 0}. The initial and boundary 
conditions for ([1]) are given as following: 

w|en\r = 0, (w^ -^70;) |r = 0, (2) 

^lan\r = 0, (^. + /3^)|r = 0, (3) 

a;|t=o = a;o(x, z), u\t=Q = UQ{yL, z) . (4) 
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This 3D model with viscosity in both u and u components was first proposed by Hon and 
Lei in [3]. The only difference between this 3D model and the reformulated Navier-Stokes 
equations is that convection term is neglected in the model. If one adds the convection 
term back to the 3D model, one would recover the full Navier-Stokes equations. This model 
preserves almost all the properties of the full 3D Navier-Stokes equations. Despite the 
striking similarity at the theoretical level between the 3D model and the Navier-Stokes 
equations, the former seems to have a very different behavior from the full Navier-Stokes 
equations. In a recent paper [5], we prove that the above 3D model with partial viscosity 
develops a finite time singularity for a class of initial condition using a mixed Dirichlet 
Robin boundary condition. 

The analysis of finite time singularity formation of the 3D model [5] uses the local 
well-posedness result of the 3D model. The local well-posedness of the 3D model can be 
proved by using a standard energy estimate and a mollifier if there is no boundary or if the 
boundary condition is a standard one, see e.g. [6]. For the mixed Dirichlet Robin boundary 
condition we consider here, the analysis is a bit more complicated since the mixed Dirichlet 
Robin condition gives rise to a growing eigenmode. 

There are two key ingredients in our local well-posedness analysis. The first one is to 
design a Picard iteration for the 3D model. The second one is to show that the mapping that 
generates the Picard iteration is a contraction mapping and the Picard iteration converges 
to a fixed point of the Picard mapping by using the Contraction Mapping Theorem. To 
establish the contraction property of the Picard mapping, we need to use the well-posedness 
property of the heat equation with the same Dirichlet Robin boundary condition as uj. The 
well-posedness analysis of the heat equation with a mixed Dirichelet Robin boundary has 
been studied in the literature. The case of 7 > is more subtle because there is a growing 
eigenmode. Nonetheless, we prove that all the essential regularity properties of the heat 
equation are still valid for the mixed Dirichlet Robin boundary condition with 7 > 0. 

2 The main result 

The local existence result of our 3D model with partial viscosity is stated in the following 
theorem. 

Theorem 2.1 Assume that uq G H^~^^{Q), loq G H'^{Q) for some s > 3/2, uo\dn = 
uozldn = and ojq satisfies (0j. Moreover, we assume that (3 G 6*00 (or S^) as defined 
in Lemma \2.1\ Then there exists a finite time T = T {\\uQ\\us+im\,\\i^o\\H'>(n)) > such 
that the system ^ with boundary condition ^,IB^ and initial data ^ has a unique solu- 
tion, -uG C{[Q,T],H'+^{n)), wGC([0,r],fl'^([7)) an(i'i/'GC([0,r],if^+2(O)). 

The local well-posedness analysis relies on the following local well-posedness of the heat 
equation and the elliptic equation with mixed Dirichlet and Robin boundary conditions. 
First, the local well-posedness of the elliptic equation with the mixed Dirichlet and Robin 
boundary condition is given by the following lemma [5]: 



Lemma 2.1 There exists a unique solution v E H''{Q) to the boundary value problem: 

-Av = f, (x,z)gJ], (5) 

v\dn\r = 0, {v, + H\r = 0, (6) 

i/ /3 G 5oo = {/3 I /3 / ^ for all k e 1?}, f G H'^'^{n) with s > 2 and f\Qn\r = 0. 
Moreover we have 

\\v\\H'>{n)<Cs\\f\\H--2(^n), (7) 

where Cg is a constant depending on s, \k\ = \/k'f + ^2- 

Definition 2.1 Let K. : H^~'^(Q,) — )■ H'^{Q,) be a linear operator defined as following: 

for all f G H'^~ ($7), /C(/) is the solution of the boundary value problem ^-(^. 

It follows from Lemma |2. II that for any / G H^^^{Q,), we have 

mf)\\HHn)<Cs\\f\\Hs-2^ny (8) 

For the heat equation with the mixed Dirichlet and Robin boundary condition, we have 
the following result. 

Lemma 2.2 There exists a unique solution oj G C{[Q,T];H^{Q)) to the initial boundary 
value problem: 

ujt = i^Au, (x, z) G fl, (9) 

wbn\r = 0, (w, +7w)|r = 0, (10) 

uj\t=o = ujo{x,z). (11) 

for lvq G H'^{Q) with s > 3/2. Moreover we have the following estimates in the case ofj > 

Mt)\\H^in) < C(7,s)e^^'*||u;o|U.(n), t > 0, (12) 

and 

Mt)\\Hsin)<C{^,s,t)\\^o\\LHn), t > 0. (13) 

Remark 2.1 We remark that the growth factor e"^^ * in U^) is absent in the case 0/7 < 
since there is no growing eigenmode in this case. 

Proof First, we prove the solution of the system ([9])- (jlip is unique. Let wi, 002^ H^{n) 
be two smooth solutions of the heat equation for < t < T satisfying the same initial 



condition and the Dirichlet Robin boundary condition. Let to = ui — lj2- We will prove 
that cj = by using an energy estimate and the Robin boundary condition at F: 



2dtJn 
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Jr 
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'dn-dz — 2uy I / uiuizdx.dz 



'dxdz -\ — / \u)z\ dyidz + 2uy I uj dxdz 
2 Jn Jn 



|V(x!| d^dz + 2vy / UJ d^dz, 



(14) 



where we have used the fact that the smooth solution of the heat equation u decays to zero 
as z — )• oo. Thus, we get 



1 d 



-— / uj^dyidz < 2uy l co^dyidz. 

Ju 



It follows from Gronwall's inequality 



g-4^7't / io'^dxdz < I u^dxdz = 0, 



(15) 



(16) 



since ojq = 0. Since oj € H^{il.) with s > 3/2, this implies that a; = for < t < T which 
proves the uniqueness of smooth solutions for the heat equation with the mixed Dirichlet 
Robin boundary condition. 

Next, we will prove the existence of the solution by constructing a solution explicitly. 
Let r]{x,z,t) be the solution of the following initial boundary value problem: 



rjt = i^Ar], (x, z) £ Q, 
V\dn = 0, r?|i=o = %(x,z), 

and let ^(x, t) be the solution of the following PDE in Ox: 

Ct = z^Ax^ + uj'^C, X G Ox, 
^l9nx = 0, ^|t=o = wo(x), 



(17) 
(18) 



(19) 
(20) 



where Ax = -4^ + -4^ and ct;o(x) = 2j f^ wo(x, z)e '^^dz. From the standard theory of the 
heat equation, we know that rj and ^ both exist globally in time. 



We are interested in the case when the initial value i]q(x, z) is related to ujq by solving 
the following ODE as a function of z with x being fixed as a parameter: 

1 _ _ 
??02 + % = wo(x, z) -a;o(x)e "^^j r?o(x, 0) = 0. (21) 

7 

Define 

uj{x,z,t) = ri, + ri + ^{x,t)e-^' , (x,z)gJ1. (22) 

7 

It is easy to check that uj satisfies the heat equation for t > and the initial condition. Ob- 
viously, uj also satisfies the boundary condition on d^\T. To verify the boundary condition 
on r, we observe by a direct calculation that {uJz + ^uj)\r = — ^??2z|r- Since r7(x, z)|r = 0, 
we obtain by using rjt = vAt] and taking the limit as z — t- 0+ that Ar/|r = 0, which implies 
that rjzzlr = 0. Therefore, u also satisfies the Dirichlet Robin boundary condition at T. 
This shows that t^; is a solution of the system ([9|)-pip. By the uniqueness result that we 
proved earlier, the solution of the heat equation must be given by (|22p . 

Since r] and ^ are solutions of the heat equation with a standard Dirichlet boundary 
condition, the classical theory of the heat equation [1] gives the following regularity esti- 
mates: 

\\v\\Hs{n] < C ||??o||H»(f^) , ||C(x)||j:^.(t^^) < Ce"''' * ||a;o(x)||^.(n^) . (23) 

Recall that ijzzlr = 0. Therefore, r]z also solves the heat equation with the same Dirichlet 
Robin boudary condition: 

{r]z)t = i^Ar]z, (x,z)gO, (24) 

{Vz)z Ir = 0, {vz)\dn\T = 0, {Vz)\t=o = Voz{x, z), (25) 

which implies that 

ll^2ll//''(n) - ^ II^Ozllfl-s(Q) . (26) 

Putting all the above estimates for r/, ijz and S, together and using (j22p . we obtain the 
following estimate: 



I'^WH^iQ) 



-Vz + V + C{^,t)e-^" 
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H^{n) 



< — \\n II -I- llnll -I- II/^i'y i-\p~'y^\\ 



It remains to bound ||%z||ji/s(q), ||%||_H-s(n) and \\^o{^)\\H''{n^) ™ terms of ||(^o||_H-»(f7)- ^y 
solving the ODE (12 ip directly, we can express rj in terms of loq explicitly 



r'OO 

r?oix, z) = -7e^^ j e^^^' f{x, z')dz' = 7 / e-^(^'-^)/(x, z')dz', (28) 



where /(x, z) = a;o(x, z) — uJo{x.)e '^^ and we have used the property that 

/>oo 

/ f{x,z)e-^'dz = 0. 
Jo 

By using integration by parts, we have 

%.(x, z) = -7/(x, z) + 7M e-^(^'-^)/(x, z')dz' = 7 / e-^(^'-^)/,,(x, z')dz'. 

J z J z 

By induction we can show that for any a = (ai, 025 03) > 

e-^("'-")L>"/(x,z')dz'. 
Let -ftr(2;) = ^e~'^^x{^) ^-iid xl-^^) be the characteristic function 



X{z) 



0, z < 0, 

1, z>0. 



Then D"r]Q can be written in the following convolution form: 

POO 

D''r]o{x,z)= K{z' -z)D''f{x,z')dz'. 

Jo 

Using Young's inequality (see e.g. page 232 of [2]), we obtain: 
WD'^voh^n) < \\K{z)\\m^+)\\D''f\\L2^n) 

Moreover, we obtain by using the Holder inequality that 

= ;^||i"°'-'-.(-.^)|L,,„,. 

Substituting dMI to ([MI) yields 

WD^mWmn) < C(7,a)(||I?"u;o|lL2(n) + ||^^"^'"^^'^o 
which implies that 

\\r]o\\H-(n) <Ci-f,s)\\ujo\\jjs^^), V s > 0. 
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1/2 



1/2 



L2{n) 



(29) 



(30) 



(31) 



(32) 



(33) 



(34) 



(35) 



(36) 



It follows from ([MD that 

ll^o(x)||j:^.(n^) < C(7)||a;o||j;^.(f^), V s > 0. (37) 

On the other hand, we obtain from the equation for tjq (|2ip that 

\\r]Oz\\H-(n)=l\\f + Vo\\H'>in)<C{-f,s)\\uo\\j^s(^), V s > 0. (38) 

Upon substituting ([36l) - (p8]) to ([27|) . we obtain 

Mln-^n) < C(7, s)e'''^'* \\^o\\H'>[n) > (39) 

where C(7, s) is a constant depending on 7 and s only. This proves (I12p . 

To prove (J13p . we use the classical regularity result for the heat equation with the 
homogeneous Dirichlet boundary condition to obtain the following estimates for i > 0: 

\\r]\\H''(n)<Cit)\\vo\\LHn), (40) 

hzWH-in) < Cis,t)\\voz\\L^(n), (41) 

M^)\\H^in^) < Cis,t)e''^'' \Poi^)\\LHn^) ^ (42) 

where C{s,t) is a constant depending on s and t. By combining (j40])-(j42]) with estimates 
(f36l) - (f38l) . we obtain for any i > that 

M\H-{n) < C{-f,s,t) (||r?02||L2(n) + WvoWmn) + e"'^^* \Po{^)\\L^n^)) 

< C7(7,s,t)||u;ollL2(c), (43) 

where C(7, s, i) < 00 is a constant depending on 7, s and t. This proves (fT3|) and completes 
the proof of the Lemma. D 

We also need the following well-known Sobolev inequality [3]. 
Lemma 2.3 Let u,v ^ H^{^) with s > 3/2. We have 

\\uv\\H'(n) <c\W\\H''{n)\\v\\H''(n)- (44) 

Now, we are ready to give the proof of Theorem 12.11 
Proof of Theorem 12.11 Let v = v?. First, using the definition of the operator K. (see 
Definition 12.1 p . we can rewrite the 3D model with partial viscosity in the following equivalent 
form: 

ivt = ^vKiu:) ^ (x,z)G = 0. X (0,00), (45) 

I Wj = Vz + u/alo 

with the initial and boundary conditions given as follows: 

w|en\r = 0, (w^ + 7^^) |r = 0, (46) 

oj\t=Q = a;o(x, z) E W, v\t=o = vo(x, z) E V'+^, (47) 



where V'+^ = {v £ H'+^ : v\dn = 0,t;^|an = 0,Vzz\dn = 0} and W = {w £ H' : w\dn\r = 
0,{wz+7w)\r = 0}. 

We note that the condition uo\oq = uoz\dn = imphes that vo\dn = ■"Ozlan = vozzlan = 
by using the relation vq = Uq. Thus we have vq G V^^^. It is easy to show by using the 
M-equation that the property uolen = uoz\dn = is preserved dynamicaUy. Thus we have 
V £ V'+^. 

Define U = {Ui,U2) = {v,uj) and X = C{[0,T];V'+^) x C7([0,T]; W) with the norm 

||C/||x = sup \\Ui\\jjs+irQ)+ sup \\U2\\jjs(Q), MU £ X 

te[0,T] tG[0,T] 

and \etS = {U£X: \\U\\x < M}. 

Now, define the map <l> : X — t- X in the following way: let ^{v,u)) = {v,uj), then for any 
t£ [0,T], 

v{x,z,t) = vo{x,z,t) +4 v{x,z,t')IC{oJ)z{x,z,t')dt', (48) 

Jo 
uj{x,z,t) = £(-u^,a;o;x,z,t), (49) 

where ti;(x, z, t) = C{vz,ujo;x, z, t) is the solution of the following equation: 

^t = Vz + z^Aw, (x, z) G ri = ilx X (0, oo), (50) 

with the initial and boundary conditions: 

'^l9Q\r = 0, {ujz + 'yuj)\r = 0, uj\t=o = ujo{x, z) . 

We use the map <l> to define a Picard iteration: U ~^^ = <!>([/) with U^ = {vq,ujq). In the 
following, we will prove that there exist T > and Af > such that 

1. U^ £ S, for all k. 

2. IIC/'^+i-C/'^IL < i||C/'=-C/^-HL, for ah k. 

Then by the contraction mapping theorem, there exists U = {v^uj) £ S such that $(C/) = U 



(51) 



which implies that C/ is a local solution of the system ()45p in X. 

First, by Duhamel's principle, we have for any g £ C([0, T]; V^) that 



£(5,a;o;x,z,t)=P(wo;0,t)+ / Vig;t',t)dt', 

Jo 



where V{g; t' , t) = g{x, z, t) is defined as the solution of the following initial boundary value 
problem at time t: 

gt = i^Ag, (x, z) G Jl = r^x X (0, oo), (52) 

with the initial and boundary conditions: 

9\dn\r = 0, {gz+ig)\r = 0, g{x,z,t') = g{x,z,t'). (53) 



We observe that g{x, z, t') also satisfies the same boundary condition as uj for any < t' < t 
since g = v^ and v'^ G y+^. 

Now we can apply Lemma 12.21 to conclude that for any t' < T and t £ [t',T] we have 



\\V{g;t',t)\\Hs^n) < C7(7,s)e'^T'(*-*')||5(x,z,t')k^(n)- 
which implies the following estimate for C: for all t G [0,^], 



(54) 



\\/:{g,u;o;yi,z,t)\\jjs(Q)<C{'y,s)e''^^^ \\\ujo\\H'{Q)+t sup \\g(yL,z,t')\\Hs^Q)] . (55) 
Further, by using Lemma [2. II and the above estimate (|55p for the sequence U 



(v^,uj^) 



we get the following estimate: 



,,fc+i 



„,+i,o^ ^ lko|lH=+i(n)+4T sup 
< Ikoll/fs+ifQ) +4r sup 

t6[0,T] 



t>'=(x,z,t) 
7;'=(x,z,t) 



sup 
sup 



IC{uj^),{x,z,t) 

Uj''{-K,Z,t) 



H^{n) 



Vt G [0, T] (56) 



u 



k+1 



< 0(^,3)6"^ * \\uJo\\Hs(n)+^ sup 
^"(f^) \ t'e[o,t] 

< C(7, s)e'''^'^ l|wo||//''(m +^ sup 



v^;{x,z,t'] 



t£[0,T] 



_H'=+i(n) 



, vtG[o,r]. 



Next, we will use mathematical induction to prove that if T satisfies the following inequality: 

8C(7,5)re-^'^(||L.o||^.(f,) + 2r||7;o|lH.+i(n)) <1 (58) 

then for all A; > and t £ [0, T], we have that 

<2\\vo\\Hs+i(n)^ (59) 



H=+'^{n) 



UJ 



wio) 



< C{j, s)e''^'"^ (\\^o\\H-{n) + 27" lbo|lH=+i(f^)) 



(60) 



First of all, U^ = {vq,ujq) satisfies ([59]) and ([60]) . Assume U^ = {v^,ui'^) has this property, 
then for f/'^+^ = (u'^+^jO;'^"'"^), using ([56]) and ([57]) . we have 



,.A:+1 



„,+i,o^ ^ lbo|lH=+i(n)+4T sup 



t>'=(x,z,i) 



sup 

^^+Hf^) tG[0,T] 



w'^(x, z,t) 



H^(Q.) 



< 2||?;o|li^.+i(f7), VtG[0,r]. (61) 



(57) 



w 



fc+i 



< 0(7,5)6^^^"^ ||a;o||H.(f^)+T sup 



/f^+l(n) 



< C(7,s)e''^'^ (||a;o||^.(f^) + 2r||z;o||//.+i(n)) > VtG[0,T] 



(62) 



Then, by induction, we prove that for any k > 0, U = {v ,io ) is bounded by (j59p and 

We want to point that there exists T > such that the inequahty (j58p is satisfied. One 
choice of T is given as following: 

8C(7, sje'"^" ( \\uJo\\Hs(n) + 2 lko|lHs+i(f7) 



Ti = min 



1 . 



(63) 



Using the choice of Tin ([63]) . we can choose M = 2 \\vo\\ jjs+i m\+C {'y , s)e'^'^ ( ||wo||j:^s(f7) + 2 ||uo|liifs+i(Q) 
then we have U^ G S, for all fc. 

Next, we will prove that $ is a contraction mapping for some small < T < Ti. 
First of all, by using Lemmas 12.11 and \'1.'6\ we have 



yk+l _ ^k 



H=+^(n) 



t j-t 



< 



+ 






H^+'^(n) 



< T sup 

te[o,T] 

+T sup 
te[o,T] 



t;'^ - /-I 



sup 

^^+Hf^) ie[o,T] 



/C(w^ 



,,fc-i 



sup 

^^+Mf^) tG[0,T] 






_f/=+l{n) 



ff^+Hf^) 



< MT sup 

\te[o,T] 



v'^ -v^ 1 



+ sup 



UJ^ - Uj'' ^ 



H'^in) 



(64) 



On the other hand. Lemma 12.21 and (j5ip imply 

£(f^,wo;x,i) - £(t!^~\wo;x,t) 



u;*^+i - 0;^= 



//'■(C) 



ii'^(n) 



< 



V{v';-vl-'-t\t)dt' 







< TC{-f,s)e''^'^ sup 

iG[0,T] 

< rC(7,s)e'^T'^ sup 

te[o,r] 






yk _yk 1 



_H'''(n) 



(65) 



Let 

T = min 



8C(7, 5)6^"^ ll^o|lH«(m + 2 ||?^o|Ih«+i 



Then, we have 



jjk+i _ jjk 



'm 



-1 r 



2C(7,s)e 



i'7 



-1 1 

' 2M' 



1 .(66) 



1 
< - 

X - 2 



U" - U' 



fc-i 
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This proves that the sequence U converges to a fixed point of the map ^ : X ^ X, and 
the limiting fixed point U = {v,uj) is a solution of the 3D model with partial viscosity. 
Moreover, by passing the limit in (j59p -(|60 p . we obtain the following a priori estimate for 
the solution v and to: 

ll^ll_H'=+i(f7) - 2 ll'"o||_H-s+i(Q) , (67) 



for < t < T with T defined in (1661) . 

It remains to show that the smooth solution of the 3D model with partial viscosity is 
unique. Let (vi,uji) and (^2,^2) be two smooth solutions of the 3D model with the same 
initial data and satisfying ||^i||j:/s+i(Q) < M and ||wj||j|^sm) < M for i = 1,2 and < t < T, 
where M is a positive constant depending on the initial data, 7, s, and T. Since s > 3/2, 
the Sobolev embedding theorem [1] implies that 

l|wj||L->(n) < \\vi\\H^+i^Q) < M, i = 1,2, (69) 

\\K^{^i)z\\L^(n) < \\^{^i)z\\H-(n) ^ ^s W^iWn^in) ^ ^sM, i = 1,2. (70) 

Let i; = f 1 — z;2 and lo = coi — u}2- Then (v, co) satisfies 

{ , / , (x,z) G Jl = J7xX (0,00), (71) 

with wlgoyr = 0, {ijJz + 7'^) |r = 0, and a;|t=o = 0, v\t=o = 0. By using ([MIl-dZQ]), and 
proceeding as the uniqueness estimate for the heat equation in (fH]) . we can derive the 
following estimate for v and w: 

J^\\v\\h^n) < C'ldl^lli^co) + M\h{n))^ (72) 

'^•' "^ <C3(|b||i.(^) + ||a;||i.(^)), (73) 



dt 
d 

^INII2^(n) 



where Ci {i = 1,2,3) are positive constants depending on M, i/, 7, Cg- In obtaining the 
estimate for (j73p . we have performed integration by parts in the estimate of the f^-term 
in the w-equation and absorbing the contribution from a;^ by the diffusion term. There is 
no contribution from the boundary term since f 1^=0 = 0. We have also used the property 
||/C(a;)2||L2(f^-) < Cs||ti;||^2(f^), which can be proved directly by following the argument in 
the Appendix of [5]. Since vq = Q and wq = 0, the Gronwall inequality implies that 
ll^llL2(n) = ||<^||l2(Q) = for < t < T. Furthermore, since v € H^~^^ and u: € H^ with 
s > 3/2, V and a; are continuous. Thus we must have v = oj = for < t < T. This proves 
the uniqueness of the smooth solution for the 3D model. D 
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